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Abstract 



We analyse the chances of detecting charged Higgs bosons of the Minimal Super- 
symmetric Standard Model (MSSM) at the Large Hadron Collider (LHC) in the 
W^h mode, followed by the dominant decay of the lightest Higgs scalar, h — > bb. 
If the actual value of is already known, this channel offers possibly the op- 
timal final state kinematics for charged Higgs discovery, thanks to the narrow 
resonances appearing around the and h masses. Besides, within the MSSM, 
the H — > W^h decay rate is significant for not too large tan (3 values, thus offer- 
ing the possibility of accessing a region of MSSM parameter space left uncovered 
by other search channels. We consider both strong (QCD) and electroweak (EW) 
'irreducible' backgrounds in the 36-tagged channel to the gg — > tbH~ production 
process that had not been taken into account in previous analyses. After a series 
of kinematic cuts, the largest of these processes is tbW^h production in the con- 
tinuum. However, for optimum tan/3, i.e., between 2 and 3, the charged Higgs 
boson signal overcomes this background and a narrow discovery region survives 
around M H ± 200 GeV. 
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1. Introduction 



The discovery of charged Higgs bosons will provide a concrete evidence of the multi- 
doublet structure of the Higgs sector. Recent efforts have focused on their relevance to 
Supersymmetry (SUSY), in particular in the MSSM, which incorporates exactly two 
Higgs doublets, yielding - after spontaneous EW symmetry breaking - five physical 
Higgs states: the neutral pseudoscalar (A), the lightest (h) and heaviest (if) neutral 
scalars and two charged ones (H^). 

In much of the parameter space preferred by SUSY, namely M#± > and 
1 < tan/3 < m t /m b p], [J, the LHC will provide the greatest opportunity for the 
discovery of particles. In fact, over the above tan/3 region, the Tevatron (Run 2) 
discovery potential is limited to charged Higgs masses smaller than m t |4j]. 

However, at the LHC, whereas the detection of light charged Higgs bosons (with 
Mh± < m t ) is rather straightforward in the decay channel t — > bH + for most tan/3 
values, thanks to the huge top-antitop production rate, the search is notoriously difficult 
for heavy masses (when M H ± > mt), because of the large reducible and irreducible 
backgrounds associated with the main decay mode H~ — > bi, following the dominant 
production channel bg — > tH~ |5|]. (Notice that the rate of the latter exceeds by far 
other possible production modes @-||, this rendering it the only viable channel at the 
CERN machine in the heavy mass region.) 

The analysis of the H — > bi signature has been the subject of many debates [fj- 
|T2"| , whose conclusion is that the LHC discovery potential is satisfactory, but only 



provided that tan/3 is small (~ 1.5) or large (~ 30) enough and the charged Higgs boson 
mass is below 600 GeV or so. 

A recent analysis [13| has shown that the rv decay mode, indeed dominant for light 
charged Higgs states and exploitable below the top threshold for any accessible tan/3 
H[ , can be used at the LHC even in the large M#± case, in order to discover H ± 
scalars in the parameter range tan (3 ~ 3 and 200 GeV < M H ± < 1 TeV. Besides, if the 
distinctive r polarisation |J| is used in this channel, the latter can provide at least as 
good a heavy signature as the H — > bi decay mode (for the large tan/? regime 

At present then, it is the tan (3 ~ 3 region of the MSSM which ought to be explored 
through other decay modes, especially those where direct mass reconstruction is pos- 
sible. The most obvious of these is the H ± — > W^*'h channel |L8| (see also ||19|), 
proceeding via the production of a charged gauge boson and the lightest Higgs scalar 
of the MSSM, with the former on- or off-shell depending on the relative values of M#± 
and Mh. In fact, its branching ratio (BR) can be rather large, competing with the 
bottom-top decay mode and overwhelming the tau-neutrino one for M H ± ~ m t at low 
tan/?: see Figs. Besides, under the assumption that the h scalar has previously 
been discovered (which we embrace here), its kinematics is rather constrained, around 
two resonant decay modes, W — > 2 jets (or lepton-neutrino) and h —>■ bb, an aspect 
which allows for a significant reduction of the QCD background. As demonstrated in 
Ref. | 20fl , signals of charged Higgs bosons in the 2 ~ tan/3 ~ 3 range can be seen in this 



channel, provided that 200 GeV ~ M#± ~ 220 GeV (see also |2T[ for an experimental 
simulation). The above lower limit on tan/3 corresponds to the border of the exclusion 
region drawn from LEP2 direct searches for the MSSM h scalar, whose mass bound is 



now set at M h £ 100 GeV or so p 



It is the purpose of this letter that of resuming the studies of Ref. HSOfl , by analysing 
the contribution to the background due to several irreducible processes, not considered 
there, whose presence could spoil the feasibility of charged Higgs searches in the W ±<y *^h 
mode of the MSSM. 

The plan of this paper is as follows. In the next Section we discuss possible signals 
and backgrounds, their implementation and list the values adopted for the various 
parameters needed for their computation. Section 3 is devoted to the presentation and 
discussion of the results. Conclusions are in Section 4. 



2. Signals and backgrounds 



We generate the signal cross sections by using the formulae of Ref. [IT]. That is, we 
implement the 2^3 matrix element (ME) for the process 

gg — > tbH~ + charge conjugate (c.c). (1) 

This nicely embeds both the gg — > tt — > tbH~ + c.c. subprocess of top-antitop produc- 
tion and decay, which is dominant for m t ~ Mjj±, as well as the bg — ► tH~ + c.c. one 
of 6t-fusion and i7 ± -bremsstrahlung, which is responsible for charged Higgs production 
in the case m t ~ M#± |23| . The ME of process ([!]) has been computed by means of the 



spinor techniques of Refs. Pi|-P7|. 

In the H W'^h -> W'^bb channel, assuming high efficiency and purity 
in selecting/rejecting 6-/non-6-jets, possible irreducible background processes are the 
following (we consider only the gg-initiated channels): 

1. the tbW~h continuum; 

2. tbW~Z production, especially when M z ~ M h ; 

3. the QCD induced case tbW~g; 

4. and, finally, tbW~H and tbW~A intermediate states; 

in which H, h, A, Z,g —>■ bb, plus their c.c. channels. Once the top quark appearing in 
the above reactions decays, two bosons are present in each event. We will eventually 
assume the W + W~ pair to decay semi-leptonically to light-quark jets, electrons/muons 
and corresponding neutrinos. Furthermore, we will require to tag exactly three b- 
jets in the final state (e.g., by using /i- vertex or high px lepton techniques). The 



same 'signature' was considered in Ref. [20|, where only the 'intrinsic' tbH — > ttbb 



background and the QCD noise due to l ti + jet' events were studied (with jet signifying 
here either a b-, light-quark or gluon jet, the latter two mistagged for the former). 
Both signal and background MEs have been integrated numerically by means of 



VEGAS [f2§] and, for test purposes, of RAMBO g9[ and Metropolis |30]] as well. While 
proceeding to the phase space integration, one also has to fold in the (x, Q 2 )-dependent 
Parton Distribution Functions (PDFs) for the two incoming gluons. These have been 
evaluated at leading-order, by means of the package MRS-LO(05A) |3l[] . 



The numerical values of the SM parameters are [I = e,/z): 



m e = m Vt = m u = m d = m s = m c = 0, 

m b = 4.25 GeV, m t = 175 GeV, 

M z = 91.2 GeV, Y z = 2.5 GeV, 

M w ± = 80.2 GeV, T w = 2.2 GeV. (2) 

As for the top width T t , we have used the LO value calculated within the MSSM (i.e., 
r t = 1.55 GeV if M H ± > m t ). 

Concerning the MSSM parameters, we proceed as follows. For a start, we assume 
that the mass of the lightest neutral Higgs particle (but not tan/?) is already known, 
thanks to its discovery at either LEP2, Tevatron (Run 2) or from early analyses at 
the LHC itself. Thus, for us, is a fixed parameter, assuming for reference the 
following discrete values: e.g., 90, 100, 110, 120 and 130 GeV. Then we express all 
other Higgs masses as a function of M H ± and tan (3. For the pseudoscalar Higgs boson 
mass, the tree-level relation Mjj± = M^ v± + M\ is assumed. Radiative corrections 
then, of arbitrary perturbative order, are in practice embedded in the H mass and 
the mixing angle a. In general, notice that, at the 'Renormalisation Group improved' 
one-loop level ||32|| , it is only for very large values of the lightest stop mass and of the 
squark mixing parameters that Mh can escape the LEP2 bound in the low tan f3 region, 
on which we will focus most of our attention. 

Finally, notice that we develop our discussion at the parton level, without consider- 
ing fragmentation and hadronisation effects. Thus, jets are identified with the partons 
from which they originate and all cuts are applied directly to the latter. In particu- 
lar, when selecting 6-jets, a vertex tagging is implied, with a finite efficiency, e&, per 
each tag. Moreover, we assume no correlations among multiple tags, nor do we in- 
clude misidentification of light-quark (including c-quark-)jets produced in W ± decays 
as 6-jets. 



3. Results and discussion 

As a preliminary exercise, we study the total production and decay cross sections before 
any cuts, as all our reactions are finite over their entire phase spaces (recall that m& ^ 0). 
This is done in Figs. §-4j for the signal and the five background processes discussed in 
the previous Section, for five values of tan/?, over the range 140 GeV ~ M H ± < 500 
GeV, for M h = 90 and 100 GeV, in the channel X -> bb, where X = h,Z,g,H or A. 
(Of course, the W ± Z and W^g backgrounds have no dependence on any of the three 
parameters aboveP].) As for the decay rates of the top (anti) quark and the boson, 
for sake of simplicity, we take them equal to 1 for the time being. The signal is always 
dominated by the QCD background and - at large tan j3 - also by the EW ones. Notice 
the local maxima of the signal rates at M#± w M\y± + Mh, as induced by the opening 
of the H~ — > W~h decay (compare to Figs. and the minima as well, due to the 

onset of the H~ —>■ bi channel instead. 



2 Note that the rates in Figs. S-0 account for the c.c. production modes as well. 



In the reminder of our analysis, we assume semi-leptonic decays of W + W pairs, as 
in Ref. [pOfl : i.e., W + W~ — > 2 jets (hereafter, jet refers to a non-fe-jet and I = e, //). 



However, as compared to that analysis, we make one simplification. Namely, we assume 
that one top (anti) quark and the IV ± boson generated in its decay have already been 
reconstructed, e.g., by using the mass selection procedure advocated in Ref. [^Cj, either 
leptonically or hadronically. This allows us to greatly reduce the complexity of our 
numerical calculation while - we believe - substantially un-affecting the relative rates 
of signal and backgrounds (in fact, all processes described produce the same final state 
and all involve at least one top quark). Then we apply the following cuts on the 
remaining particles (here, the label j refers to the decay products of the second W ± 
boson present in the event, which can be either light-quarks or leptons): 

Pt(&J, missing) > 20 GeV (3) 

on the transverse momentum (including the missing one), 

[77(6,j)| < 2.5 (4) 

on the pseudorapidity, and 

Ai?(66,6j,jj) > 0.4 (5) 
on the relative separation of h- and light-quark jets/leptons j, where 



AR(ij) = ^Arj(ij) 2 + A0(*7) 2 , (6) 

is defined in terms of relative differences in pseudorapidity T](ij) and azimuth 4>(ij), 
with i j = b, ]/$.. Furthermore, we impose (see also Ref. [pOj] ) 



\M bb — Mh\ < 10 GeV (7) 

on exactly one pair of fe-jets, 

Mjj > 50 GeV (8) 

on the light-jet (or lepton- neutrino) pair (recall that the W can be off-shell), and, 
finally, 

\M bbm -M t \ <20 GeV (9) 

around the top mass if three 6's are present in the event (in addition to the one already 
used to reconstruct the top (anti)quark). In such a case, one may assume that the 
charged Higgs boson has predominantly been produced in the decay of a top (anti) quark 
(when M H ± ~ m t ). If instead only two appear, then one should conclude that the Higgs 
has mainly been generated in a bremsstrahlung/fusion process (because M#± ~ m t ) 
with a 6-(anti)quark lost along the beam pipe. Our 2^3 production mechanism 
naturally allows one to emulate both dynamics in a gauge invariant fashion, including all 
interference effects. As already mentioned, however, we will assume a triple 6-tagging, 
this implying an overall efficiency factor of e b multiplying our signal and background 
rates. (Thus, the third 6-jet in eq. (|]) is actually non-6-tagged: it can be interpreted as 
the jet system satisfying neither eq. (0) nor eq.(H).) We take e b = 0.5, like in |2y, plf| 
(and assume 100% lepton identification efficiency). 



Given the signal production rates before acceptance and selection cuts, it is clear 
that - for such an e& - even at high collider luminosity (i.e., / Ldt = 100 fb -1 per 
annum), hopes of disentangling the charged Higgs boson of the MSSM in the W^*'h 
decay channel are only confined to the very low tan/3 region. We will thus restrict 
ourselves to study in the reminder of the paper tan (3 values which are, e.g., below seven. 
The total signal rates after the cuts have been applied can be found in Fig. |^, 

for the choices tan/3 = 1, 2, 3 and 7, as a function of M H ±. For reference, we illustrate 
the 'borderline' case M h = 100 GeV. (Indeed, a lower M h value at tan /3 = 2 is in 
contradiction with LEP2 data, whereas higher masses induce a far too large suppression 
on BR(/f ± — ► W^h): see Fig. ^.) The trends in the figure are the consequence of two 
effects. On the one hand, the production cross section of gg — > tbH~ + c.c. is roughly 
proportional to (mf cot /3 2 +m^ tan /3 2 ), so that its maxima occur at very low or very high 
tan/3. On the other hand, we have seen how the largest H — > W~^h decay fraction 
is attained for tan/3 ~ 2. In the end, the largest values for a(gg — > tbH~) x BR(H~ — > 
W~^*'K) + c.c. are obtained for tan/3 = 1: see Fig. [5|. Unfortunately, such a tan/3 
value is already excluded in the MSSM from LEP2 data [20]. For the optimal remaining 
choice, i.e., tan/3 = 2, the annual rate never exceeds 140 events (before any 6-tagging 
efficiency but after acceptance and selection cuts). The maximum occurs at M H ± w 200 
GeV, significantly above the real threshold at Mh + M^± pa 180 GeV. 

We now compare such a signal with the irreducible backgrounds 1.-4., for the same 
choice of tan/? and Mh (where relevant). This is done in the upper half of Fig. [!], at 
the level of total production rates. After the cuts ©-((J) are enforced, all background 
components in 2.-4. are overwhelmed by the signal in the vicinity of Mjj± = 200 GeV, 
whereas the W ± h continuum production is always larger than the H — > W h resonant 
channel. Thus, it is relevant to compare the last two processes in the 'reconstructed' 
invariant mass M w ± h , i.e., that obtained from pairing the two 6-jets fulfilling condition 
(0) and the two light-quark jets (or, alternatively, the lepton-neutrino pair) satisfying 
eq. (§) and not already reconstructing M w ± on their own and m t in association with any 
of the 6's (see Ref. [pOfl ). The spectrum in this variable is presented in the lower half of 
Fig. §, for our ideal case M H ± = 200 GeV (and, again, tan/3 = 2 and M h = 100 GeV). 
For such MSSM parameter combination, the charged Higgs signal is well above the 
continuum for values of M^± h which are ±20 GeV from Mh±. (To vary M^± and/or 
tan p basically corresponds to rescale the solid line in the last plot by a constant factor, 
according to the rates in Fig. |5|.) 

For reference, Tab. [I] presents the number of events of resonant and continuum W ± h 
production at the LHC, after 300 inverse femtobarns of collected luminosity, for e\ = 
0.125, in the window \M H ±-M w ± h \ < 40 GeV, for the three values M H ± = 180, 200 and 
220 GeV. Although very small, a H — > W~h signal is generally observable above the 
W~h continuum for Mjj± around 200 GeV. Our numbers are roughly consistent with 
those in Ref. [pC[, if one considers that we neglect the finite efficiency of reconstructing 
one boson and the associated top (anti) quark and since we have chosen somewhat 
different cuts. Therefore, in the end, the dominant backgrounds remain (in the 36- 
tagged channel) the H~bi + c.c. decay and the QCD noise involving misidentified 



gluons, i.e., those already identified in Ref. [20 



4. Conclusions 



In summary, in this paper, we have complemented a previous analysis of the pro- 
duction and decay of charged Higgs bosons of the MSSM at the LHC, in the channels 
gg — > tbH~ and H — > W~(*'h (and charged conjugated modes), respectively, by 
considering several irreducible backgrounds in the 36-tagged channel, i.e., tbH~ + c.c. 
— > 3b 2 jets £ + 'missing energy' (where the initial 6-(anti)quark is usually lost along 
the beam pipe), which had not yet been considered. 

We have found that, after standard acceptance cuts and a kinematic selection along 
the lines of the one outlined in Ref . p(| , the dominant background among those con- 
sidered here is the continuum production gg — > tbW~^h + c.c. However, the latter 
has been found to lie significantly below the signal in the only region where this is de- 
tectable: when tan/5 ^2 — 3 and M H ± rs 200 GeV (with M h around 100 GeV, close to 
the latest LEP2 constraints). Thus, the chances of detecting the H~ — > W~h — > W~bb 
decay in such a (narrow) region of the MSSM parameter space depend mainly on the 
interplay between this mode, the competing one H — > bt — > W~bb and the QCD 
background with mistagged gluons, as are the latter two that clearly overwhelm the 
former (recall the last figure in ||2Q|| ). 

We have carried out our analysis at parton level, without showering and hadronisa- 
tion effects but emulating typical detector smearing. We are confident that its salient 
features should survive a more sophisticated simulation, such as the one presented in 
Ref. |2"I|| . Besides, our results concerning the backgrounds can be transposed to the 
case of non-minimal SUSY models (where the H ± discovery potential can extend to a 
much larger portion of parameter space), such as those considered in Ref. p0| , so that 
also in these scenarios the irreducible backgrounds analysed here can be brought under 
control. 
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Table 1: Number of signal (S, tbH~) and dominant background (B, tbW~h) events, 
along with the statistical significance S/ \/B, after the implementation of the cuts (|3|)- 
(P) in the decay channel W + W~h — > 2 jets I i>a bb. Rates are given for tan/3 = 2, 
M h = 100 GeV, three choices of M H ±, as obtained by using the MRS-LO(05A) set of 
PDFs, after 300 fb -1 of luminosity and for = 0.5. 
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Figure 1: Dominant branching ratios of the charged Higgs boson of the MSSM for 
selected values of tan (3 over the mass range 140 GeV ~ M#± ~ 500 GeV. The mass of 
the lightest Higgs boson has been fixed at Mh = 90 GeV. 




Figure 2: Branching ratios into W h pairs of the charged Higgs boson of the MSSM 
for selected values of tan /3 over the mass range 140 GeV ~ M#± ~ 500 GeV. The mass 
of the lightest Higgs boson has been fixed at M h = 90, 100, 110, 120 and 130 GeV. 
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Figure 3: Cross sections of signal and backgrounds for selected values of tan (3 over the 
mass range 140 GeV ~ M H ± £ 500 GeV. The mass of the lightest Higgs boson has been 
fixed at M h = 90 GeV. Here, both the top quark and the W ± boson are kept on-shell 
and no decay rates and cuts are applied. 




Figure 4: Same as Fig. g for M h = 100 GeV. 




Figure 5: Cross sections of signal and backgrounds for selected values of tan (5 over the 
mass range 160 GeV ~ M H ± ~ 500 GeV, after the cuts (|f)-(|9]) and including decay 
rates. The mass of the lightest Higgs boson has been fixed at M h = 100 GeV. 
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Figure 6: (Top) Total cross sections of signal and backgrounds for tan (3 = 2 and Mh = 
100 GeV after the cuts and including decay rates, as a function of the charged 

Higgs boson mass over the range 160 GeV ~ M H ± ~ 500 GeV. (Bottom) Differential cross 
sections in the reconstructed W^h invariant mass, for a M^± = 200 GeV signal and for 
the dominant background after the cuts ©-© and including decay rates. Gaussian 
smearing of all transverse momenta is included: with {&{pt)/pt) 2 = (0.60/^/pr) 2 + 
(0.04) 2 for jets and (<j(pt)/pt) 2 = (0.12/ \fp~r) 2 + (0.01) 2 for leptons/missing particles. 



